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One-particle interchain hopping in a system of coupled Luttinger liquids is investigated by use of 
exact diagonalizations techniques. Firstly, the two chain problem of spinless fermions is studied in 
order to see the behaviour of the band splitting as a function of the exponent a which characterizes 
the ID Luttinger liquid. We give numerical evidence that inter-chain coherent hopping (defined 
by a non-vanishing splitting) can be totally suppressed for a ~ 0.4 or even smaller a values. The 
transverse conductivity is shown to exhibit a strong incoherent part even when coherent inter-chain 
hopping is believed to occur (at small a values). Implications for the optical experiments in quasi-lD 
organic or high-T c superconductors is outlined. 



Recently, the study of strongly correlated fermions con- 
fined to coupled chains has received a great deal of in- 
terest in particular as a way of studying the dimensional 
cross-over from ID Luttinger-like behaviour to 2D. 

Some time ago, Anderson emphasized the crucial dif- 
ference between in-plane and inter-plane (c-axis) trans- 
port observed in copper oxide superconductors [0]. In- 
deed, experimentally the transverse conductivity has a 
completely incoherent frequency dependence S: there 
seems to be no sizeable Drude-like term (except in the 
optimally doped systems) and <j{u>) is a very slowly in- 
creasing function of the frequency. This phenomenon 
has been interpreted as an incoherent hopping or as the 
"confinement" of the electrons inside the weakly coupled 
planes 0. However, for coupled Fermi liquids (FL), Lan- 
dau theory predicts coherent transverse hopping and no 
anomalous transport. Therefore, these data have sug- 
gested that the ground state (GS) of the two-dimensional 
(2D) plane itself is not of the usual FL type. Unfortu- 
nately, it has been impossible yet to prove the NFL na- 
ture in 2D except for unrealistic models. 

However, it is well known that the generic features of 
correlated ID electrons are not FL-like but rather those of 
a Luttinger liquid (LL) Q and the precise nature (asymp- 
totic behaviour, exponents, etc..) of the system can be 
easily controlled. In addition, quasi-lD systems are re- 
alized in nature, and the problem of coupled chains is of 
direct relevance there. For instance, in the case of the 
organic superconductors of the (TMTSF^X family ||, 
also known as the Bechgaard salts, the high temperature 
properties are believed to be essentially one-dimensional, 
while the low-temperature behaviour is rather two- or 
three-dimensional. This cross-over is presumably respon- 
sible for the anomalies observed in the temperature de- 
pendence of several quantities (such as the 2k p contri- 
bution to the relaxation rate M , the ratio of the perpen- 
dicular conductivity to the parallel one ||, the plasma 
edge when the electric field is polarized perpendicular to 
the chains ra],. . . ) as well as for the insulating behaviour 
reported for (TMTSF^CIO^ in the presence of a strong 
enough magnetic field ||. Theoretically a lot of work 
has already been devoted to that problem but several as- 



pects of this cross-over have to be understood better, in 
particular those concerned with the transport properties 
perpendicular to the chains. 

Hence, from now on and for sake of simplicity, we shall 
only deal with weakly coupled chains. 

The effect of single-particle transverse hopping has pre- 
viously been studied from a renormalization group point 
of view H . Let us recall here that a LL has a different 
structure from FL: there are no quasi-particle like exci- 
tations but instead collective modes (charge and spin) 
with different velocities which lead to the so-called spin- 
charge separation; moreover, the density of states n{k) 
has no step-like structure at the Fermi level but instead 
a power law singularity n(k) — n(kp) ~ \k — kp\ a defining 
the parameter a which depends on the intra-chain inter- 
action. It turns out that the hopping t± is a relevant 
perturbation when a < 1 ||. 

However, it has been argued that relevance in that 
sense was not necessarily a sufficient condition to cause 
coherent motion between chains. This, e.g., can be seen 
from the following model fTcj ]; let a system of two sep- 
arated chains be prepared at time t = with a differ- 
ence of A7V particles between the two chains. Then, the 
interchain hopping is turned on and one considers the 
probability of the system of remaining in its initial state, 
P{t). Coherence or incoherence can then be denned as 
the presence or absence of oscillations in P(t). 

In Ref. pi , the authors found two regimes for P(t) 
which depend entirely on the value of a: the case a < 
1/2 exhibited coherent motion while a > 1/2 showed no 
signal of coherence (see also discussion in |llj). 

This paper is devoted to the study of various aspects 
of interchain coherence in systems of strongly correlated 
spinless fermions. A more extended version of this work 
can be found in Ref. 0]. We shall derive several quan- 
tities sensitive to the coherent /incoherent nature of the 
hopping transverse to the chains from exact diagonaliza- 
tions of small systems by the Lanczos algorithm jl2| . 

We consider here a model of spinless fermions on a 
lattice formed by m chains of length L with a weak inter- 
chain hopping: 
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where (3 labels the chain (f3 = l,...,m), j is a rung 
index (J — 1, . . . , L), Cj t p is the fcrmionic operator, and 
V(5) is a repulsive interaction between two fermions at 
a distance S (the lattice spacing has been set to one). 
For convenience, we choose a repulsive interaction of the 
form V(i) — 2V/(i + l) for i < i , with, more specifically, 
io = 1, 2, 3 which corresponds to an interaction extending 
up to first, second and third nearest neighbours (NN) 
respectively. 

In the x and y-directions, we shall use arbitrary bound- 
ary conditions (BC) by threading the system with a mag- 
netic flux & x and & y respectively (except for m — 1 or 2 
chains where open BC are used in the y-direction) which 
is measured in unit of the flux quantum <E> = he/ e. 

The motivation to introduce flux is two-fold. Firstly, 
as proposed by Kohn |l4| transport properties of a corre- 
lated system can be directly measured from the response 
of the system to a twist in the boundary condition. Sec- 
ondly, our ultimate goal is to extract quantities in the 
thermodynamic limit from finite size scaling analysis. It 
turns out that a simple way to improve the accuracy for 
a fixed system size is to average over the boundary con- 
ditions e.g. over <& y and/or <& x JlllJlrj] . 

In order to characterize the behaviour of coupled 
chains, it is required to, first, compute the parameters 
of a single chain for the same model. Indeed, one impor- 
tant issue is to study whether interchain transport is a 
universal function of the LL parameters only or whether 
it depends on the details of the model. In the case of NN 
interactions, the hamiltonian (known as the t-V model) 
can be mapped onto a spin chain problem by a Jordan- 
Wigner transformation and this is exactly soluble by the 
Bethe ansatz; thus, a is known for each filling ||. How- 
ever, for extended interactions in space, a numerical in- 
vestigation is necessary with the help of conformal in- 
variance identities. It turns out that the exponent a can 
be related to simple physical quantities |l7| which can 
be easily extracted from standard exact diagonalization 
results using the Lanczos algorithm. 

The simplest approach to investigate interchain coher- 
ence is to consider two coupled chains, i.e. a 2 x L ladder. 
We proceed along the lines of Ref. (ll| . In the absence of 
interaction, t± leads to bonding and anti-bonding dis- 
persion bands corresponding to transverse momentum 
k± = or 7r respectively, as seen in Fig. ([j]). The splitting 
2t± between these bands can be viewed as the signature 
of a coherent transverse hopping. These bands corre- 
spond to a (5-function singularity in the single particle 
hole (electron) spectral function for k < kp (k > fop). 

In the case of interacting particles, this (5-function sin- 
gularity is replaced by a power law singularity and the 
elementary excitations are collective modes. Here, we 
address the issue of the influence of the hopping t±_ on 



this singularity, in particular we investigate whether a 
splitting occurs. 




FIG. 1. Dispersion relation along the chain direction in the 
open shell configurations. Full (open) symbols correspond to 
occupied (empty) states and \x is the chemical potential. 

As can be seen in Fig. ([!]), in the open shell configu- 
ration, one can add or remove a particle exactly at the 
Fermi momentum which is independent of the system 
size. We thus define the splitting as the difference be- 
tween the electron and hole excitation energies i.e. 

AE = E (N e + l,k e ) +E Q (N e - l.kfc) - 2E Q {N e ). (1) 

where Eo(N e ) is the reference energy corresponding to 
the absolute GS of the N e — nN electron system. The 
momenta for the electron and hole excitations are fixed, 
k e = (kp,7r) and k/j = (kp,0). For V=0, this expres- 
sion exactly gives the splitting 2t± for any system size. 
For t± = but finite interaction strength, AE is finite, 
however it scales to zero in the thermodynamic limit as 
expected. We observe too that an accurate finite size 
scaling analysis can be performed for finite interaction 
strength and finite t±, assuming 1/L finite size correc- 
tions. 

The extrapolated values (L = oo) of AE are plotted as 
a function of t 1 /^ ^ as suggested by the renormalization 



group approach [ 
in Fig. (ft. 



for various interactions (io = 2 model) 
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FIG. 2. Extrapolated (L = oo) values of AE vs t 1 /^ q) 
in the io = 2 spinless fermion model with several values of V 
(indicated on the figure). 
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A strong reduction of this splitting for increasing a in- 
dicates that intrachain repulsion has a drastic influence 
in prohibiting interchain coherent hopping. For large in- 
teractions like V = 3.25 (a ~ 0.38), there is a critical 
value t*^ (a) of t± below which incoherent transverse hop- 
ping takes place. This corresponds to the case a > olq. 
For smaller interaction like V = 1 (a ~ 0.08) our data 
are consistent with a coherent behaviour when t± — > 0. 
Hence, transverse hopping remains coherent in this case 
which corresponds to a < a.$. For intermediate interac- 
tions like V — 2 (a ~ 0.22), our data are not conclusive. 
However, very small values of «o like 0.2 (or even smaller) 
are not inconsistent with our numerical analysis. 

When spin is taken into account, the spin-charge sepa- 
ration that occurs in ID should suppress even further the 
coherent transverse hopping. At a qualitative level, this 
can be understood from the fact that only real electrons 
can hop from one chain to the next and this is believed 
to become more difficult in the presence of spin-charge 
separation Although dealing with a different filling 
n = 1/3 and with particles with spin, the same qualita- 
tive behaviour is found numerically in Ref. ]l3[ |; however, 
the decrease of the splitting with a is stronger. 

The previous study suggests that the interaction tends 
to confine the electrons within the chains, although no 
complete confinement seems to occur at small a values 
where AE ^ 0. A better understanding of this phe- 
nomenon can be achieved by investigating the transport 
properties along the y-axis (inter chain) and more pre- 
cisely the transverse optical conductivity which is the 
linear response of the system to a spatially uniform, time 
dependent electric field in the transverse direction. For 
such a study, a torus geometry is needed (to > 3) so that 
a current can flow around the loop in the y-direction. 
One of the main advantage of the optical conductivity is 
that it can directly be measured experimentally 

The real part of the optical conductivity can be written 
as a sum of two parts, a yy (uS) — 2irD yy S(u>) 
where 2irD yy is the Drude weight in the y-direction 

Note that in order to check the numerical results, we 
can use an important f-sum rule [fl9| : the integrated con- 
ductivity is proportional to the mean value of the trans- 
verse kinetic energy in the ground state. 

As originally noted by Kohn jl4j , D yy can be obtained 
from the dependence of the ground state energy Eq on 
$ y as 
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total kinetic energy remains finite down to vanishing t± 
as in the non-interacting case. We first consider a fixed 
value of the flux $ y such that ||p = 0. In that case, fi- 
nite size effects are found to be already weak for the two 
largest cases L — 8 and L = 12 and the Drude weight is 
strongly suppressed compared to the V = case [Q . 
In order to mimic the case of many parallel chains, 



an average over $ y is realized 
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and I„ 



do 



i.e. we calculate 
which is obtained 



by averaging the transverse kinetic energy. We observe 
that the values obtained for the Drude weight by aver- 
aging the 3x8 data over $ y are very close to the ones 
obtained on the 3 x 12 cluster at constant flux. Note 
that, however, if the number of coupled chains is kept 
fixed (here m = 3), even in the limit L — > oo, we ex- 
pect (D vy )$ 7^ D yy {<& y ). In any case, when confinement 
within each of the individual chains starts to occur we do 
not expect crucial differences between the cases of 3 or of 
an infinite number of coupled chains (if only the t± term 
couples the chains). In Fig. |3|, we plot the Drude weight 
and the total sum rule as a function of t± for a moder- 
ate interaction. The behaviour of these two quantities is 
not incompatible with the t\ law of the non-interacting 
case |2(]]. However, the intrachain interaction has drastic 
effects. Firstly, the total sum rule is strongly reduced 
compared to V = (shown as a reference on Fig. ||). 
Secondly, it is found that n{D yy )^ y and I yy behave dif- 
ferently. 
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where N = mL is the number of sites. 

The previous quantities have been calculated numeri- 
cally on finite 3 x L lattices using the Lanczos algorithm. 
We have chosen a quarter filled band so that the extrap- 
olation of results for L = 4, 8 and 12 is possible. 

By choosing adequate boundary conditions along x 
open shells can be realized as seen in Fig. (|l|). In this 
case, even for a finite system, the Drude weight and the 



FIG. 3. K{D yy )<s> y (■) and I yy (•) as a function of t± for 
a 3 x 8 system and V — 2, io = 3. For comparison, the V = 
case is also displayed as a dotted line. 

Let us consider the ratio r — ■n(D yy )^, y / I yy which cor- 
responds to the relative part of the Drude weight in the 
total optical conductivity. In Fig. ^, we note that the 
ratio r decreases as t± goes to and can become rather 
small. 
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FIG. 4. Ratio of the Drude weight to the total conductivity 
in the 3x8 system as a function of tx for NN and longer range 
interaction with the same NN magnitudes V = 2. 

These results should have very important consequences 
on the experimental side. Indeed, our results predict, 
for small t±, anomalous transport perpendicular to the 
chains even when D yy does not completely vanish (if 
a < «q) since spectral weight is suppressed predomi- 
nantly from the coherent part of the conductivity. In or- 
der to confirm this behaviour, the frequency dependent 
optical conductivity has been calculated directly by use 
of the Kubo formula, 
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where j y is the transverse current operator and the sum 
runs over all the excited states. 
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FIG. 5. Transverse optical conductivity vs frequency for 
a 3 x 8 system at quarter filling with V — 2, 2q = 3 and 
t± = 0.1. The Drude 5-function has been represented with 
the same small imaginary part e = 0.04. 



In the free case, the current operator commutes with 
the Hamiltonian and therefore, the conductivity only 
contains a Drude peak. However, as can be seen on Fig. ||, 
for finite interaction strengths and in the 3x8 cluster, a 
pronounced structure appears at finite frequency. In or- 
der to determine more accurately the position in energy 
of the weight, we have computed the first moment of the 
distribution, 
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which is expected to behave smoothly with the various 
parameters. We observe that, once t± is turned on be- 
tween the chains, weight immediately appears predom- 
inantly at finite frequencies (see Fig. |^). This typical 
frequency increases with the strength of the interaction. 
This is clearly a signature of some form of incoherent 
perpendicular transport. 
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FIG. 6. First moment of the conductivity as a function of 
t± for V = 2 and different interaction ranges. 

In this work, different approaches to interchain coher- 
ence have been investigated. As a first step, we have fo- 
cused on the energy splitting generated by the transverse 
hopping in the dispersion relation of the LL collective 
modes. By finite size scaling analysis, we have shown 
that this splitting was monitored by the LL parameter 
a. However, incoherent interchain hopping is found for 
much smaller values of a than those predicted by the RG 
calculations Moreover, in the regime where tj_ is still 
relevant (a < ao), the most important results are that 

(i) the Drude weight and the total optical sum rule grow 
less rapidly with t± than in the non interacting case and 

(ii) even when the Drude weight remains finite (when t± 
is relevant), transverse transport is predominantly inco- 
herent in the small tj_ regime. How small tj_ needs to be 
so that this regime is observed depends on the strength 
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of the interaction. Typically, for a ~ 0.2, strong sup- 
pression of coherent transport occurs up to t± ~ 0.15. 
This phenomenon could explain the anomalous transport 
which is observed experimentally. 
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